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ABSTRACT 

A  model  of  resonance  fluorescence  of  a  laser-driven  two-level 
atom  near  a  metal  surface  is  reviewed,  and  results  are  shown  in  the 
weak-field  limit  for  the  population  inversion  and  the  power  spec¬ 
trum.  Extension  of  the  model  to  interacting  adatoms  is  discussed 
in  terms  of  phonon-induced  level  broadening  and  adatom  migration. 


INTRODUCTION 

Recently  there  has  been  considerable  theoretical  and  experi¬ 
mental  interest  in  the  interaction  of  electromagnetic  radiation 
with  surfaces.  The  perfection  of  the  fatty-acid  monolayer  assembly 
technique  has  led  to  a  series  of  experiments  in  which  the  fluores¬ 
cence  of  an  excited  atom  or  molecule  at  a  fixed  distance  from  a 
metal  surface  (gold,  silver  and  cooper)  was  measured. I  Experi¬ 
ments  have  also  demonstrated  enhanced  luminescence  and  resonance 
fluorescence  for  adspecies  on  metals  (and  also  on  dielectric 
fibers)  as  compared  to  the  surface-free  molecules .2 *3  in  order 
to  understand  the  excitation-relaxation  behavior  and  resonance 
fluorescence  spectrum  of  a  laser-driven  atom  near  a  metal  surface, 
we  choose  in  this  paper  a  simple  two-level  model  for  the  adatom, 
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combined  with  the  techniques  of  quantum  optics.  We  shall  first 
review  the  model  for  a  single  adatom  and  then  extend  it  to  many 
interacting  adatoms. 


SINGLE  TWO-LEVEL  ADATOM 


We  assume  the  adatom  to  be  driven  by  a  partially-coherent 
laser  field.  Although  the  atom  has  no  dipole  moment  in  its 
ground  state,  it  can  have  a  significant  induced  transition  dipole 
in  a  driving  field.  The  reflected  field  provided  by  a  metallic 
mirror  and  surface  plasmon  resonance  can  influence  the  dynamic 
behavior  and  light  scattering  spectrum  of  the  adatom.  The  field 
reflected  by  the  metallic  surface  takes  the  quantum  form4 

ER(t)  =  y12C12{t)  f(d)  +  y21a21(t)  f*(d)*  (1) 

where  y  is  the  transition  dipole  between  adatomic  states  |l>  and 
|j>,  o  *J  =  |  i><j |  is  the  adatomic  transition  operator,  and  f(d)  is 
a  frequency-dependent  and  adatom-metal  distance-dependent  (d)  func¬ 
tion  which  will  be  discussed  in  more  detail  later. 


In  the  rotating-wave  approximation  (RWA) ,  the  surface-dressed 
optical  Bloch  equations  (SBE)  can  be  written  (assuming  the 
equilibrium  condition  W  *  -1) 


r  ' 

u 

f  4 

CM 

CO 

-Y2  +  iA  iff  (t)/2  0 

s2i(t) 

W(t) 

= 

ifi+(t)  -Yi  -i^'(t) 

W(t) 

0  -ifi+(t)/2  -Y2-iA 

s12(t) 

~  y  i 


(2) 


W(t)  is  the  population  inversion;  S^jCt)  =  (t)  exp  ( -iuc  (t) ;  = 

n  exp[+i4>)t)  ] ,  where  =  2 1 y„,  |E0(t;/fi  is  the  regular  Rabi  frequency 
and  the  electric  field  E(t)  of  laser  is  written  in  terms  of  the 
amplitude  E  (t)  and  phase  factor  $  as  E  (t)  exp  (-ia^t  -  i<t>(tj  +  c.c.; 
A  =  ui2^“w^°is  the  detuning,  where  is  the  transition  frequency 
of  the  adatom  and  u,  is  the  laser  frequency;  the  surface-induced 


^L  _ 

phase-decay  constant  is  yg  =  2  Im(f)  | y.  _  |Vfi;  y^  =  Yo  +  Y  »  Y1  and 
are  the  population  andphase  decay  constant,  respectively ,  in 
the  surface-free  case.  We  have  implicitly  included  the  interatomic 
collision  effects  by  letting  y^  t  2y2» 
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To  calculate  yg,  we  must  know  the  complex  function  f(d),  which 
can  be  determined  by  the  Sommerfeld-Hertz  vector  procedure, ^ 


Im(f) 


1+E 


1  ,  2  w 

(w  -  U)  )  + 

sp 


,2  4 

0  ID 

_ E£_ 


4,2 

oj  6 


{■—■  +  ~^~2~  ~  sinD - ^  cosD] 


2  1 

2cj6uj  [n  cosD  +  — -r  sinD]  }  > 
P  D 


(3) 


where  £.  =  1  and  n  =  — -  +  —  for  the  case  of  the  induced  dipole 
oriented  parallel  toDthe  surface,  £  =  2  and  n  =  for  the  per¬ 
pendicular  case,  e^  is  the  dielectric  constant  for  the  gas  medium, 
kj  *=  is  the  wave  number,  to  is  the  emission  frequency,  and 

D  =  2kjd  is  the  reduced  distance  which  is  dimensionless.  The 
Drude  model  £2(1*3)  “  1  -  — ^as  ^een  utilized,  where  6 
is  the  inverse  of  the  relaxation  time,  <0p  is  the  plasma  frequency, 
and  io_p  *  10  (1+ej)  l'2  is  the  surface  plasmon  frequency  of  the 
interface. 


Equation  (3)  leads  to  a  surface-induced  phase-decay  constant 
Yg  (in  the  unit  of  Einstein's  spontaneous  decay  constant  A)  of  the 
following  form: 


?a  fa. (n 
o  1 


sinD  - 


cosD 


)  + 


a2(h 


cosD  + 


sinD 


)], 


(4) 


where  a0  =  3/8[(l-B2)2  +  (a/<oSp)  ],  a,  *  (2-B2)2  +  (2a/(osp)  -  B4, 

&2  m  (2B>2a/iog  ,  and  B  =i^/^sp,  a  =  6/6  and  ej  -  1  has  been  assumed. 
The  above  exprlssion  for  y8  can  be  further  simplified  by  assuming 
perfect  reflection,  i.e.,  R4  ■  R11  ■  -1,  where  R"*-  (K**)  is  the  reflec¬ 
tion  coefficient  for  an  incident  ray  polarized  perpendicular  (paral¬ 
lel)  to  the  plane  of  incidence.  Then 


3  „  ,  .  n  cosD.  , 

Y  ■  rC  (n  sinD - 5^) .  (5) 

D  6 

We  use  the  phase-diffusion  model  for  the  laser  field: 
«fi"(t1)J2+(t2)»  *  fi2  exp(-  1 1 1 1),  where  is  the  laser  band¬ 
width.  We  consider  the  weak-field  or  large-detuning  limit,  wherein 
W(t)  »  <<W(t)>>  «  -1.  The  double  bracket  signifies  an  average  over 
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the  stochastic  ensemble  and  a  quantum  mechanical  average.  Assuming 
the  laser  field  to  be  turned  on  at  t  =  0  with  a  constant  amplitude 
for  t  >0,  we  can  obtain  an  analytical  form  for  the  population 
inversion. ^  The  power  spectrum  S(w)  as  a  function  of  the  emission 
frequency  u>  is  given  by  the  Fourier  transform  of  the  dipole-dipole 
correlation  function,  and  for  steady  state  it  can  also  be  obtained 
analytically  to  yield  the  expression^ 


S(u))  - 


(y2+yl)2+a2 


{• 


(y2-yl)2+a2 


l 


YL^Y2~YL+^  +  2yL^ 


Y?  + 


*?2^2~yL+^2)  +  2''LA(a,_ul21)  2^Y2+YL^Y2 

2  2  ^  +  2  2~ 
Y2  +  (u)-u21>  yi  Iy2  +  (u)-<d21)  ] 


(6) 


Sample  results  for  these  analytical  forms  of  W(t)  and  S(w)  are 
shown  in  Figure  1.  The  time  evolution  of  W(t)  is  shown  in  the 
upper  part  of  the  figure  for  two  different  adatom-surface  plasmon 
resonance  conditions:  g  *  0.5  and  g  *  0.9.  We  have  used  the  value 
of  3.5  x  10"3  for  the  ratio  d/wSp»  which  corresponds  to  a  silver 
surface.  We  make  the  following  three  observations:  (i)  Each  curve 
approaches  a  stationary-state  value  for  large  t.  (ii)  The  value  of 
W  depends  on  the  adatom-surface  plasmon  resonance  condition  g  and 
the  reduced  adatom-surface  distance  D.  Generally,  the  values  of  W 
for  the  off-resonance  case  (g  =  0.5)  are  less  than  those  for  the 
near-resonance  case  (g  =  0.9).  Furthermore,  the  values  of  W  for 
larger  D  are  less  than  those  for  smaller  D,  since  the  levels  of  the 
adatom  for  smaller  D  are  strongly  broadened  by  the  adatom-surface 
interaction,  so  that  the  population  of  the  upper  level  increases, 
(iii)  For  large  D  (curve  3),  the  curves  for  g  «*  0.5  and  g  *  0.9 
coincide,  since  the  adatom-surface  coupling  tends  to  vanish.  The 
strong  oscillations  in  this  nearly  surface-free  situation  are 
caused  by  beating  between  the  driving  laser  field  and  the  "adatomic 
oscillator". 

The  bottom  part  of  Figure  1  displays  the  resonance  fluorescence 
power  spectrum,  with  the  Rayleigh  peak  near  u>  «  and  fluorescence 
peak  near  u  =  The  latter  peak  is  strongly  dependent  on  D  and 

g.  As  D  decreases  and  g  approaches  unity,  the  upper  adatomic  level 
broadens,  and  hence  the  power  in  the  fluorescence  peak  increases. 


MANY  ADATOMS 

So  far  we  have  studied  the  SBE  and  resonance  fluorescence  for 
a  single-adatom  system  in  which  the  main  effects  of  the  surface  are 
characterized  by  the  reflected  field  and  the  adatom-plasmon  inter- 
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action.  For  a  system  consisting  of  many  adatoms  as  depicted  in 
Figure  2,  the  SBE  may  be  derived  by  taking  an  ensemble  average 
over  the  adatoms.  We  note  that  the  cooperative  feature  of  a  many- 
adatom  system  can  be  significantly  different  from  that  of  the 
single-adatom  system  due  to  static  processes,  such  as  dipole-dipole 
coupling,  and  dynamical  processes,  such  as  collisional  broadening 
among  the  mobile  adatoms.  We  also  note  that  the  adatoms  together 
with  the  surface  atoms  may  actually  form  a  "pseudomolecule"  whose 
vibrational  degrees  of  freedom  may  be  coupled  to  the  electronic 
excitations.  To  describe  the  many-body  features  of  the  system,  one 
can  write  a  total  Hamiltonian  as 


H  = 


ha  +  hb  +  hab  +  haa  +  haba  +  har  +  hbb  +  haf* 


(7) 


and  Hg  are  the  unperturbed  Hamiltonians  for  the  adatoms  (A) 
and  the  phonon  bath  modes  (B) ,  respectively,  with  a  coupling 
interaction  H  H  and  H  represent  the  adatom-adatom  inter¬ 
actions  via  direct  clipole-aipole  coupling  or  phonon-mediated 
coupling;  Har  represents  the  interaction  between  the  adatoms  and 
reflected  field  and/or  surface  plasmon;  Hgg  describes  energy  dif¬ 
fusion  on  and  into  the  substrate  via  phonon-phonon  coupling;  and 


LASER  RADIATION 


Figure  2.  Schematic  diagram  for  a  many-adatom  system  in  which  two- 
level  adatoms  are  excited  by  laser  radiation  and  coupled  to  the 
surface  plasmon  and  substrate  phonons  as  described  by  the  total 
Hamiltonian  of  Equation  (7)  . 
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H  _  Is  the  radiative  coupling  between  the  applied  laser  field  and 
the  adatoms.  To  derive  the  SBE,  a  second-quantized  form  of  the 
above  Hamiltoniam,  involving  dipole  operators  for  the  adatoms  and 
harmonic  operators  for  the  phonon  modes,  is  appropriate.  The  mathe¬ 
matical  form  may  be  found  in  References  7  and  8  and  is  not  shown 
here.  Instead,  we  want  to  discuss  below  the  physical  mechanisms 
associated  with  the  individual  terms  of  the  total  Hamiltonian. 


From  the  SBE  we  see  that  the  overall  level  broadening  of  an 
adatom  is  characterized  by  phase  relaxation  (T^)  and  energy  relaxa¬ 
tion  (Fi).  The  relaxation  factors  may  be  decomposed  as 


ri  =  Yi  +  Y1B  +  Y1r 


P2  =  y2  +  Y2B  +  y2R  ‘ 


(8) 

(9) 


Y j  and  Y2  are  the  previously-defined  relaxation  factors  or  decay 
constants  for  a  single-adatom  system  except  that  the  collision- 
induced  broadening  in  is  modified  further  by  the  migration  of 
the  adatoms.  Y^g  and  Y2B  are  the  relaxation  factors  due  to  electron- 
phonon  coupling,  and  y^r and  y^d  are  due  to  the  reflected  field. 
Depending  on  the  interaction  mechanisms,  these  phase  and  energy 
relaxation  factors  are  influenced  differently  by  the  system's  para¬ 
meters  such  as  the  adatom-surface  distance,  the  coverage  or  the 
adatom-adatom  distance,  the  surface  temperature  and  the  frequency 
spectra  of  the  whole  system.  Some  of  the  important  features  are 
the  following: 


(i)  The  phonon-induced  broadening  Yj>  resulting  from  H  and 
is  strongly  temperature  dependent  through  the  phonon  occupation 
number  and  is  governed  by  the  frequency  spectra  of  the  phonon  modes. 
For  example,  for  an  adatom  with  frequency  closer  to  the  Debye  fre¬ 
quency,  we  expect  a  stronger  coupling,  or  larger  yb* 

(ii)  The  dephasing  broadening  ^  due  to  the  adatom-plasmon 

interaction  may  be  further  broadened  via  and  in  which  the 

collisional  effects  may  be  significant  for  a  high-coverage  system 
and/or  highly  mobile  adatoms  In  physisorbed  states. 

(iii)  The  ref lected-f ield  broadening  y^r  is  related  to  Re(f) 
and  Y2R  t0  Io(f)  due  to  H^.  The  direct  (H.^)  and  phonon-mediated 
(H.-.T  adatom-adatom  interactions  through  tneir  near  dipole  field 
ana  exchange  of  photons  are  slowly-varying  functions  of  their 
separations  compared  with  the  dependence  of  those  effects  due  to 

H._  on  the  reduced  distance. 

AK 

We  finally  note  that  the  influence  of  the  surface  on  the  SBE 
and  the  resonance  fluorescence  is  exhibited  through  the  level 
broadenings  characterized  by  the  above  relaxation  factors  Tj  ^  and 
the  adatomic  frequency  shift.  To  include  the  effects  of  the’ 
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frequency  shift,  the  detuning  A  defined  earlier  should  be  replaced 
be  an  effective  detuning  A.  The  effective  resonance  condition  A  = 

0  may  be  met  by  means  of  the  continuum  phonon  spectra.  Furthermore, 
when  a  two-level  adatom  moves  close  to  the  surface  to  form  a  bond, 
e.g.,  physisorbed  to  the  substrate,  it  becomes  a  multilevel  system 
in  which  each  electronic  level  has  many  vibrational  sublevels. 


SUMMARY 

A  model  has  been  developed  to  describe  the  dynamical  behavior 
of  a  two-level  adatom  at  a  metal  surface.  Surface-dressed  optical 
Bloch  equations  (SBE),  which  account  for  surface-reflected  photons 
and  adatom-surface  plasmon  coupling,  have  been  derived.  Collisions 
with  foreign  gas  atoms  and  effects  of  the  laser  bandwidth  have  been 
included.  The  population  inversion  and  the  power  spectrum  of  the 
scattered  light  for  a  single  adatom  have  been  evaluated.  These 
quantities  exhibit  behavior  which  is  strongly  dependent  on  the 
adatom  surface  distance  and  adatom-surface  plasmon  resonance.  For 
a  many-adatom  system,  the  SBE  and  the  resonance  fluorescence  are 
influenced  not  only  by  static  interactions,  such  as  dipole-dipole 
coupling,  electron-phonon  coupling  and  the  phonon-mediated  adatom- 
adatom  interactions,  but  also  by  dynamical  type  processes,  such  as 
laser-enhanced  migration  of  the  adatoms. 
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